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Abstract 

Single pseudo-scalar and vector mesons inclusive photoproduction 7/1 — > 
MX via higher twist mechanism is calculated using the QCD running cou- 
pling constant method. It is proved that in the context of this method a higher 
twist contribution to the photoproduction cross section cannot be normalized 
in terms of the meson electromagnetic form factor. The structure of infrared 
renormalon singularities of the higher twist subprocess cross section and the 
resummed expression (the Borel sum) for it are found. Comparisons are made 
with earlier results, as well as with leading twist cross section. Phenomeno- 
logical effects of studied contributions for w, K, p-meson photoproduction 
are discussed. 
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1 INTRODUCTION 



One of the fundamental achievements of QCD is the prediction of asymptotic scaling 
laws for large-angle exclusive processes and their calculation in the framework of 
perturbative QCD (pQCD) [1-3]. In the context of the factorized QCD an expression 
for an amplitude of an exclusive process can be written as integral over x, y of 
hadron wave functions (w.f.)Q $j(x, Q 2 ) (an initial hadron), $*(x, Q 2 ) (a final 
hadron) and amplitude T H (x, y; as(Q 2 ), Q 2 ) of the hard-scattering subprocess [2]. 
The hard-scattering amplitude T# (x, y; as(Q 2 ), Q 2 ) depends on a process and can 
be obtained in the framework of pQCD, whereas the w.f. $(x, Q 2 ) describes all 
the non-perturbative and process-independent effects of hadronic binding. The 
hadron w.f. gives the amplitude for finding partons (quarks, gluons) carrying the 
longitudinal fractional momenta x = (xi, x 2 , ■■■X n ) and virtualness up to Q 2 within 
the hadron and, in general, includes all Fock states with quantum numbers of the 
hadron. But only the lowest Fock state (qiq 2 - for mesons, uud - for proton, etc.) 
contributes to the leading scaling behavior, other Fock states' contributions are 
suppressed by powers of l/Q 2 . In our work we shall restrict ourselves by considering 
the lowest Fock state for a meson. Then, x = x%, x 2 and x% + x 2 = 1. 

This approach can be applied for investigation, not only exclusive processes 
but also for the calculation of higher twist (HT) corrections to some inclusive pro- 
cesses, such as large-pT dilepton production [4], two-jet +meson production in the 
electron-positron annihilation [5], etc. The HT corrections to a single meson inclu- 
sive photoproduction and jet photoproduction cross sections were studied by various 
authors [6,7]. In these early papers for calculation of integrals over x = X\, x 2 , like 

/~ / a 5 (Q 2 )$(x, Q 2 )F(x, a s (Q 2 ), Q 2 )5{1 - Xl - x 2 )d Xl dx 2 (1) 

which appear in an expression of the amplitude, the frozen coupling constant ap- 
proximation was used. Some comments are in order concerning this point. It is 
well known [8], that in pQCD calculations the argument of the running coupling 
constant (or the renormalization and factorization scale) Q 2 should be taken equal 
to the square of the momentum transfer of a hard gluon in a corresponding Feyn- 
man diagram. But defined in this way, as{Q 2 ) suffers from infrared singularities. 
Indeed, in a meson form factor calculations , for example, Q 2 equals to x\y\Q 2 or 
x 2 y 2 Q 2 , — Q 2 being the four momentum square of the virtual photon. In the single 
meson photoproduction jh — > MX, this scale has to be chosen equal to — X\u or 
x 2 s, where u, s are the subprocess's Mandelstam invariants [6]. Therefore, in the 
soft regions x\ — > 0, y\ — > 0; x 2 — >• 0, y 2 — ► or x\ — > 0, x 2 — > integrals (1) diverge 

1 Strictly speaking, $jvf( x :Q 2 ) is a hadron distribution amplitude and it differs from a hadron 
wave function; the former can be obtained by integrating the corresponding wave function over 
partons' transverse momenta up to the factorization scale Q 2 . But in this paper we use these two 
terms on the same footing. 
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and for their calculation some regularization methods of as(Q 2 ) in these regions 
are needed. In the frozen coupling approximation these difficulties were avoided 
simply by equating Q 2 to some fixed quantity characterizing the process. In form 
factor calculations this is Q 2 = Q 2 /4 [9], in the single meson photoproduction - 
Q 2 = S/2,-u/2 [6]. 

Recently, in our papers [10,11] devoted to the investigation of the light mesons 
electromagnetic form factors, for their calculation we applied the running coupling 
constant method, where these singularities had been regularized by means of the 
principal value prescription [12]. In our recent work we consider the inclusive pho- 
toproduction of single pseudoscalar and vector mesons 7/1 — > MX using the same 
approach. 

2 CALCULATION OF THE HIGHER TWIST 
DIAGRAMS 

The two HT subprocesses, namely r yq 1 — > Mq 2 and 7g 2 — ► Mq x contribute to the 
photoproduction of the single meson M in the reaction 7/1 — > MX . The Feynman 
diagrams for the first subprocess are shown in Fig.l. We do not provide the set of 
diagrams corresponding to the second subprocess ^q 2 - ► Mq^, they can be obtained 
from Fig.l by exchanging the quark and antiquark lines. The momenta and charges 
of the particles in question are indicated in Fig. 1(a). In our investigation the 
meson mass is neglected. As is seen from Fig.l, in the HT subprocess the meson 
M is coupled directly to the photon and the hadron quark and its suppression in 
comparison with leading twist subprocesses is caused by a hard gluon exchange in 
the higher twist diagrams. 

The amplitude for the subprocess 751 — > Mq 2 can be found by means of the 
Brodsky-Lepage method [2], 

M = [ [ dx 1 dx 2 5(l -xi - x 2 )T H (xi,x 2 ;as(Q 2 ), s, u, t)§ M (xi,x 2 ;Q 2 )(2) 
Jo Jo 

In (2), Th is the sum of graphs contributing to the hard-scattering part of the 
subprocess, which for the subprocess under consideration is 7 + qi — > (qiq 2 ) + q 2 , 
where a quark and antiquark from the meson form a color singlet state (q{q 2 ). 

The important ingredient of our study is the choice of the meson model w.f. $m- 
In this work we calculate the photoproduction of the pseudoscalar (pion, kaon) 
and vector (p-meson) mesons. For these mesons in the literature [3], [13] various 
w.f. were proposed. Here for the pion and kaon we use the phenomeno logical 
w.f. obtained in [3] by applying the QCD sum rules method, for calculation of 
the p-meson photoproduction we utilize both w.f. found in [3] and derived in [13]. 
The reason is that, in accordance with results of Ref.[13], the wave functions of 
longitudinally and transversely polarized p-mesons are similar (coincide in shape), 
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whereas in [3] a significant difference between them was predicted. In [13] the 
authors suggested also that the change in shape of the transverse p-meson w.f. may 
increase the rate of the production of transversely polarized p-mesons by a factor 
2. We think that the single meson photoproduction is a suitable arena for checking 
this conclusion. 

The pion and p-meson wave functions have the form 

^ M (x^l) = ^ sy (x)[a + b(2x-ir}. (3) 

For the model w.f. the coefficients a, b take the following values: 
Chernyak-Zhitnitsky w.f. [3]; 

a = 0,6 = 5, for the pion, 

a = 0.7, b = 1.5, for the longitudinally polarized pi — meson, (4) 
a = 1.25, b = —1.25, for the transversely polarized pr — meson. 

Ball-Braun w.f. [13]; 

a = 0.7, b = 1.5, 

for both longitudinally and transversely polarized p-meson. Here we have denoted 
by x = Xi the longitudinal fractional momentum carrying by the quark within the 
meson. Then, x 2 = 1 — x and X\ — x 2 — 2x — 1. 

The pion and p-meson w.f. are symmetric under replacement X\ — x 2 <-> x 2 — x±. 
But the kaon w.f. is non-symmetric; $k(xi — x 2 ) ^ ^k(x 2 — xi) [3]. Indeed, the 
kaon w.f. includes a term proportional to odd power of (2x — 1), 

*k(x,i$) = ^ sy (x)[a + b(2x-l) 2 + c(2x-lf], (5) 
a = 0.4, 6 = 3, c=1.25, 

and may be written as the sum of the symmetric <& s (x,/Aq) and antisymmetric 
® a (x,pl) parts, 

$,(x, pi) = ®« y (x) [a + b(2x - If] , $ a (x, nl) = <5>« y (x)c(2x - l) 3 . (6) 
In (3), (5), (6) $%,y( x ) * s ^ ne asymptotic w.f. 

0) = ^Mi-4 (7) 

where fu is the meson decay constant; f w = 0.093 GeV, fx = 0.112 GeV. In the 
case of the p-meson we take f^ = fj = 0.2 GeV for the CZ w.f., and f^ = 0.2 GeV, 
fT = 0.16 GeV for BB w.f. 

The normalization of $m(^,A*o) a ^ A*o — 0.5 GeV is given by the condition 
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1 -^,pl) = 0=. 
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The factor \[2 appearing in the normalization of a vector meson is included in the 
p- meson decay constant. 

The formalism for calculation of the HT subprocess cross section is well known 
and described in [6,14]. We omit details of our calculations and write down the 
final expression for da HT /dt. We find: 

for the pseudoscalar and longitudinally polarized vector mesons, 



da HT (e 1 ,e 2 ) 32ir 2 C F a E \ e\ 



dt 



9s 2 



s 2 



Ilt-2h(hs + hu)l + ll^ 



K\t- 2K X {K x u + K 2 s) i + K 2 L 



(9) 



2eie 2 
sut 



hKj 2 - h(K 2 s + Kttys - K^hs '+ I 2 u 



)«]}. 



for the transversely polarized vector meson, 
da HT (e u e 2 ) Un 2 C F a E -i 



dt 



9s 4 u 2 



[eiul 2 - e 2 sK 2 ] 



(10) 



In (9), (10), cue — 1/137 is the fine structure constant, C F = 4/3 is the color 
factor. The Mandelstam invariants for the subprocess are defined as 



s = (zp + q) 2 = zs, 
t = (q-P) 2 =t, 

\2 



(11) 



U 



{zp — P) = zu, 



where s, t, u are the Mandelstam invariants for the process 7/1 — * MX , z is the 
longitudinal fractional momentum of the quark qi out of the hadron h. 

The main problem in our investigation is the calculation of quantities I± t2 , K^ 2) 



h = 



1 ri 



Jo 



ff 

Jo Jo 



and 



K 2 



Jo 



Jo 



dxidx 2 5{\ - 


- Xi - 


- x 2 )a 5 (Q 2 )<l) M (a; 1 , x 2 ; Q 2 ) 


dx\dx 2 5{l - 


- Xi - 


x 2 

- x 2 )a s {Q 2 1 )^M{x 1 ,x 2 ; Q\) 






x x x 2 


1 dx\dx 2 8{l 


- x 1 


- x 2 )a s (Ql)^M(x 1 ,x 2 ; Ql) 


1 dx\dx 2 8{l 


- x 1 


Xi 

- x 2 )a s (Q 2 2 )® M (xi,x 2 ;Q 2 2 ) 



x 2 



(12) 
(13) 

(14) 
(15) 
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where for Ji, I 2 the renormalization and factorization scale is Q\ = x 2 s, for K\, K 2 
it is given by Q 2 = —x%u. 

Let us first consider the frozen coupling constant approximation. In this ap- 
proximation we put Q\ 2 equal to their mean values s/2, —u/2 and remove as(Q\ 2 ) 
as the constant factor in (12-15). After such manipulation, the integrals (12-15) are 
trivial and can easily be found. For the mesons with symmetric w.f. we get 

*© ■ 2 '*(fH*(fMf 
<-l) - »?(-D -<*(-!)'«(-! 

where superscript " 0" indicates that the quantities I, K are found in the frozen 
coupling approximation. Here the function Im(Q 2 ) is 



/m(Q 2 



1 dx$ M (x,Q 2 
x(l — x) 



In this approximation, using the last expressions and (9), (10), (12-15) one can easily 
reproduce results of [6] for the subprocess cross section^. 
In the case of the kaon we find 



1 2 



2a f 



dx§ s (x, s/2) f 1 dx§ a (x, s/2) 



1 — x 
1 dx$ s (x, s/2) 





+ 



x 



(16) 



1 — x 



It is evident that 



I° 2 (s/2) + 2J?(S/2). 



The same is also true for K\ and K\. This means that in the case of a pseu- 
doscalar meson with the non-symmetric w.f. the result of [6] is not valid and in the 
calculations our expressions (9), (12-15) have to be applied. 

The important problem in the single meson inclusive photoproduction is the 
possibility of normalization of the HT subprocess cross section (9), (10) in terms of 
the electromagnetic form factor Fm{Q 2 ) of the corresponding meson. 

The electromagnetic form factor Fm{Q 2 ) of the meson M is given by the ex- 
pression 



F M {Q 2 



Jo 



$* M (y, Q 2 )T } H f {x, y- a s (Q 2 ), Q 2 )$ M (x, Q 2 )dxdy. 



(17) 



2 The difference between our expressions and corresponding formulas in Ref.[6] is caused by our 
definition of the antiquark's charge, i.e. in our expressions the charge of the antiquark from M is 
— e2, whereas in Ref.[6] it is denoted by ei. 
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Here 



Ti f (x,y;a s (Q 2 ),Q 2 ) 



167TC; 



a s (Q 2 (l-x)(l-y)) a s (Q 2 xy) 
e i r-, rr~, \ e 2' 



(1 -*)(!-!/) 



For the meson with symmetric w.f. using the frozen coupling approximation (Q 2 
Q 2 /4) we get 



f m {Q 2 ) 



lQuC F a s {Q 2 /A) 

Q 2 



ei 



■ es) u — r3 



(18) 



It is not difficult to conclude that for such mesons the quantities (I ) 2 and (X ) 2 
in the cross section can be expressed in terms of F M 



J?(?/2) 



3(3/2) 



16ttC f 

2 a s (s/2) 



4vrCV 



[Q 2 \f m (Q 2 



Q 2 F M (Q 2 ) 



(19) 



Q 2 =2s 



&=2s 



For mesons with non-symmetric w.f. from (17) we find 



F M {Q 2 ) 



167rCW(Q 2 /4) 

Q 2 



ei - e 2 



$ s (a;,Q 2 /4)c^ 



+ 



x 



ei - e 2 J 



//■i $ a (x,g 2 /4)rfa; V 

2(ei+e2) U — r^^J ' 



l-x 



(20) 



It is now clear that (J ) 2 , (if ) 2 (16) are not proportional to F M (20). This means 
that even in the context of the frozen coupling approximation the HT subprocess 
cross section may be normalized in terms of the meson form factor only if the 
photoproduction of the meson with symmetric w.f. is considered. 
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3 THE RUNNING COUPLING CONSTANT METHOD 
AND IR RENORMALONS 



In this section we shall calculate the integrals (12-15) using the running coupling 
constant method and also discuss the problem of normalization of the higher twist 
process cross section in terms of the meson electromagnetic form factor obtained in 
the context of the same approach. 

As is seen from (12-15), in general, one has to take into account not only the 
dependence of a.{Q\ 2 ) on the sca le Q121 but also an evolution of <&m{%, Q\ 2) with 
Q\ 2 . The meson w.f. evolves in accordance with a Bethe-Salpeter type equation, 

but its dependence on Q 2 is mild and may be neglected by replacing $m(^, Q12) ~ > 
^m(x,/Mq). Such approximation does not change considerably numerical results, 
but phenomenon considering in this article (effect of infrared renormalons) becomes 
transparent. 

Let us clarify our method by calculating the integral (12); the quantities J 2 , K 12 
can be worked out in the same way. For the mesons with symmetric w.f. Eq.(12) 
in the framework of the running coupling approach takes the form 

= f 1 «s((l - x)s)$ M {x, nl)dx ^ 
J 1 

The ag((l — x)s) has the infrared singularity at x — > 1 and as a result integral (21) 
diverges (the pole associated with the denominator of the integrand is fictitious, 
because <Ev ~ (1 — x), and therefore, the singularity of the integrand at x — 1 is 
caused only by a s ((l —x)s)). For the regularization of the integral let us relate the 
running coupling at scaling variable as((l — x)s) with the aid of the renormalization 
group equation in terms of the fixed one as(s). The renormalization group equation 
for the running coupling a(s) = as(s)/ir 

^ = 1W* (22) 
has the solution [12] 

a(Xs) ~ . r - (23) 

In (22), (23), the one-loop QCD coupling constant 0:5 (/i 2 ) is defined as 

Po — 11 — 2n//3 being the QCD beta-function first coefficient. 
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Having inserted (23) into (21) we get 

I 1 (s) = a s (s) - w y , , tt, 24 

KJ byj Jo (l-x)(l + (l/i)ln(l-x))' V ; 

where t = 4n/asis)Po- 

The integral (24) is, of course, still divergent, but now it is recasted into a form, 
which is suitable for calculation. Using the method described in details in our work 
[10] it may be found as a perturbative series in a. sis) 



/if^EK^ Sn, S n = C n ®T\ (25) 



The coefficients C n of this series demonstrate factorial growth C n ~ (n — 1)!, which 
might indicate an infrared renormalon nature of divergences in the integral (24) and 
corresponding series (25). The procedure for dealing with such ill-defined series is 
well known; one has to perform the Borel transform of the series [15] 



u n-i 



£['i]W = £(^T)! C - ^ 

then invert B[Ii\iu) to obtain the resummed expression (the Borel sum) for his). 
This method is straightforward but tedious. Therefore, it is convenient to apply 
the second method proposed also in our work [11], which allows us to bypass all 
these intermediate steps and find directly the resummed expression for his). For 
these purposes let us introduce the inverse Laplace transform of 1/ (t + z) 

1 r°° 

= / exp[-(* + z)u]du. (27) 

t + z Jo 

Then h(s) may be readily carried out by the change of the variable x to z = ln(l— x) 
and using (27) 



his) 



4V37T/, 



M 



A) 

86 

i — u 



cxp 



Attu 



a s is)(3 



a + 5b 



1 — u 2 — u 



+ 



46 



u 



(28) 



Eq.(28) is nothing more than the Borel sum of the perturbative series (25) and the 
corresponding Borel transform is 



„ rr . . , a + b a + 56 86 46 
B[h]{u) = T^l- 2^7 + 3^-4^? 



(29) 
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The series (25) can be recovered by means of the following formula 



d 



n— 1 



C »=U B|/il( " ; 



«=0 



The Borel transform B[Ii](u) has poles on the real u axis at u = 1; 2; 3; 4, which 
confirms our conclusion concerning the infrared renormalon nature of divergences 
in (25). To remove them from Eq.(28) some regularization methods have to be 
applied. In this article we adopt the principal value prescription [12]. We obtain 



[h (s)Y 



M 



A) 

L*(A 3 ) 



. Li(X) . rl ,Li{\ 2 ) 
[a + 6) — ^ - (a + 56) ;„ J + 



A 



A 2 



86 



A 3 



46 



Li(\ 4 ) 
A 4 



(30) 



where Li(\) is the logarithmic integral [16], for A > 1 defined in its principal value 

(31) 



Li(\) = P.V. r^-, A = s/A 2 
Jo mx 

For other integrals from (13-15) we find 



[h (s) 



M 



[a + 6) — ^ - 46- 



A 



A 2 



46 



L?(A 3 
A 3 



and 



\k x {-u)r = ih 



-u 



, [K 2 (-u)} res = [I 2 



-u 



(32) 



(33) 



From (30), (32), (33), we conclude that in the framework of the running coupling 
approximation even for mesons with symmetric w.f. we have 



[h (ST* ^ [ h (s)} 7 



[K 2 (-u)} res * [K, {-u)\ 



Therefore only our results for the subprocess cross section (9), (10) are correct. 

Another question is, as we have discussed in Sect. 2, the normalization of the 
meson photoproduction cross section in terms of the meson elm form factor. The 
pion and kaon form factors have been calculated by means of the running coupling 
approach in our previous papers [10, 11]. Let us write down the pion form factor 
obtained using the pion's simplest w.f., that is, the asymptotic one (a = 1, 6 = in 
(4)) 



Q 2 F W (Q 2 ) 



res _ (lQirf w y 

asy (3 



- + (In A - 2) -U- + (In A + 2) 



• (34) 
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From (30), (34) it follows that the relations (19) do no longer hold. The same is also 
true for the pion's other w.f., as well as for p L — and px~ mesons. In other words, in 
the running coupling approach the HT subprocess cross section (9), (10) cannot be 
normalized in terms of the meson form factor neither for mesons with symmetric 
w.f. nor for non-symmetric ones. 

Let us, for completeness, write down I(s),K(— u) calculated for non-symmetric 
w.f. (5) 



[h(s)r = 



M 



A) 



(a + b + c) — ^—^ — ( a 



X 

Lt(X 4 ) 



Li(A 5 



2(46 + 9c)^-4( & + 5c)^ + 8c 



A 2 

Lt(X 5 ) 



+ 



A 5 



(35) 



ih (?)r 



4V37T/ 



M 



A) 



, Li(X) 0/0J n , Li(\ ) 
[a + b + c) — ^ - 2 (26 + 3c) — ^+ 
A A z 



4(6 + 3c)— 8c- 



X 3 A 4 



(36) 



The expressions for [Ki (— u)] res and [K 2 (—u)] res may be obtained from (35), (36) 
by c — > — c, A = s/A 2 — > —u/A 2 replacements, respectively. With these explicit 
expressions and the results of [11] at hand one can check our statements concerning 
the normalization of the subprocess cross section for kaons. 

Some comments are in order concerning these results. First of all, it is instructive 
to compare sources of the infrared renormalons in our case and in other QCD 
processes considered in [17]. In these articles the running coupling constant method 
was used in one-loop order calculations for resummation of any number of fermion 
bubble insertions in the gluon propagator. This technique corresponds to partial 
resummation of the perturbative series for a quantity under consideration. Indeed, 
for such quantities, the coefficients S n in (25) have the following expansion in powers 
of fa, 

S n = C n f3ft + Rnfio + . . . . 

Therefore, by defining the Borel transform as in (26) and inverting it one obtains 
a partially resummed expression for a physical quantity. In our calculations of 
the HT cross section we use the leading order term for T H ; there are no gluon 
loops in the corresponding Feynman diagrams in Fig.l. It is not accidental that 
S n (25) in our case has exactly ~ (3q~ 1 dependence and, hence, the expressions 
(30), (32), (35), (36) are exact sums of the corresponding perturbative series. In [17] 
it has been demonstrated that the only source of terms of order ~ ccg/^o" 1 in (25) 
is the running coupling constant. The source of these terms in the HT calculations 
is also the running coupling ots(s(l — x)) (or as(—u(l — x))), which runs due to 
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the integration in (1) over the meson quark's (antiquark's) longitudinal fractional 
momentum, but not because of a loop integration. 

Another question commonly discussed in papers involving IR renormalons is an 
ambiguity produced by the principal value prescription used for the regularization 
of divergent integrals (28). The ambiguity introduced by our treatment of (28) is 
a higher twist and behaves as A 2 /Q 2 (the first renormalon pole is u — 1). But the 
subprocess under consideration itself is already the higher twist one. Therefore, we 
can safely ignore such "HT-to-HT" corrections. 

At the end of this section let us write down the HT correction to the single 
meson photoproduction cross section by taking into account both HT subprocesses; 
751 — » Mq 2 and 7^2 Mq\. It is not difficult to prove that the second subprocess 
cross section can be obtained from (9), (10) by e\ <-> e 2 replacement. Then the HT 
correction to the single meson photoproduction cross section is given by 

91,52 I 

gt(*-,-* r r"^ — • (37) 



dp^dy 




where 



z = 



p T e y 



y/s - p T ey ' 

Here the sum runs over the hadron's quark q 1 and antiquark q 2 flavors. In (37) 
q^(z*,—t), q 2 (z*,— t) are the quark and antiquark distribution functions, respec- 
tively. All r.h.s. quantities are expressed in terms of the process cm. energy y/s, 
the meson transverse momentum px and rapidity y using the following expressions 



1= -, t = -pWse , u — -= -. (38) 



Eq.(37) is the final result which will be used later in our numerical calculations. 

4 PHOTOPRODUCTION OF MESONS AT THE 
LEADING TWIST LEVEL 

In our study of the single meson photoproduction a crucial point is the comparison 
of our results with leading twist (LT) ones. This will enable us to find such domains 
in the phase space in which the higher twist photoproduction mechanism is actually 
observable. 

The LT subprocesses, which contribute to a meson photoproduction are: 
a photon-quark (antiquark) scattering 

7(9) + 9i(Pi) -»■ QM + gipz), {lli -> Qi9), (39) 
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and photon-gluon fusion reactions 



l(q)+9(pi)^qi(p2) + qi(p3). (40) 

In this article we consider the inclusive cross section difference in the photon-proton 
collision, namely 



A « = - M+X) - MTy^ - M ' X) s E « + ' S «- ' < 41) 



The LT subprocess which dominates in this difference is jq — > gq with q — > M. Its 
cross section at the tree level is well known, 



dcr LT 8ira E eg 



t - s 

(*s{s)t + a s (-t)j 
s t 



(42) 



dt 3s 2 
where the Mandelstam invariants of the subprocess are 

s = (q + pi) 2 , i=(q-p 2 ) 2 , U = (q-P3) 2 - 

In (42), the running coupling constant 0:5 is evaluated at momentum scales s and 
\t\, which are equal to off-shell momenta carried by the virtual quark propagators 
in the corresponding Feynman diagrams of the leading twist subprocess •yq — * gq. 

Many other subprocesses contribute to the meson photoproduction, among them 
i) 7^ — > gq with g — > M, ii) yq — > gq with q — > M, iii)7<? — > gq with g — > M, 
lY )l9 ~^ QQ with q — > M or g — > M. Considering the cross section difference 
we not only reduce the number of subprocesses contributing to A M , because the 
subprocesses involving gluons or antiquarks contribute equally to and S^l 
and cancel in Ajff, but also solve two other important problems. The first one 
is the next-to-leading order correction to the meson photoproduction cross section 
calculated in [18]. In this paper the authors have investigated the ratio 

C7T+ _ w - = da HO ( 1P -> n+Xydprdy - da HO { 1P -> ^X)/d VT dy =R _ 1 
da Born (jp -> 7T+X) /dp T cfa/ - da Born (yp -> n-Jf) /dp T cfo/ 

for the cross section difference as a function of at ^/s = 14.1 GeV and y = 0.5. 
They have found that this ratio is negative and almost constant with p T {p T = 
2 — 6 GeV/c). This means that the K-factor for the cross section difference is 
less than 1. In other words, using the LT cross section (42) for calculations of 
Rm =| A^/A^f I a t the same or slightly different kinematic regimes, we only 
underestimate the ratio Rm and related quantities and give lower bounds for them. 
The second problem solved by our choice of A M is a contribution to the photo- 
production cross section originating from the photon's quark and gluon content. 
It is well known that the photoproduction process yp — > h + X may proceed via 
two distinct mechanisms; photon can interact either directly with the hadron's 
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partons (direct photoproduction), or via its quark and gluon content (resolved pho- 
toproduction). As was demonstrated in [18], the contribution from the resolved 
photoproduction almost completely cancel in the n + — n~ difference. These results 
obtained in [18] for pions at certain kinematic domain seemingly are valid also for 
other light mesons at the same or slightly different kinematic conditions. 

Then the leading twist contribution to the single meson photoproduction in 
7P — > MX is given by the expression, 

da LT _ ^ r 1 dxq p (x, -i)D M/q (z, -i) da LT ^ 



dPxdy g Jxmin z dt 

where 

Xy/S V s V s _ Pt£ V 

In (43), q p (x,—i) and DM/ q {z,—t) are a quark q distribution and fragmentation 
functions, respectively. The subprocess invariants s, t, u in (43) are functions of 

s, pt, y, 

f PTy/se~ y A xp Ty /se y 
s = xs, t = , u = . (45) 

z z 

Eq.(43) together with (37) for the HT contributions will be applied in the next 
section for numerical calculations. 



5 NUMERICAL RESULTS 

In this section we compute the ^p — > M + X and ^p — > M~X inclusive cross sections 
E M +, as well as the difference A M = S M + — E M - by taking into account the 

dominant LT (75 — > gq with q — > M), and HT (75 — * Mq) contributions to the 
inclusive photoproduction. Only the HT cross section of K~ photoproduction is 
calculated using the proton s and u quarks induced subprocesses, which contribute 
at the same order. Our calculations are performed for M = n , K, p at y/s = 
14.1 GeV, 25 GeV. 

In this work, for quark distribution functions, we borrow the leading order 
parametrization of Owens [19]. This parametrization is suitable for our purposes, 
because the HT mechanism probe the quark distribution functions at 
z* = Pt exp (— y) / — p^ exp (y), which for chosen process's parameters is always 
more than 0.01. The same is also true for x min in the LT cross section (43). That 
is, kinematical conditions allow us to avoid the region of small x < 0.01, where 
Owen's parametrization may give incorrect results. The same reason will enable us 
to compute S^" ignoring a contribution from the leading twist subprocess 75 — > qq, 
which otherwise may be considerable. 
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The quark fragmentation functions are taken from [20]. Recently, in [21], a new 
set of fragmentation functions for charged pions and kaons, both at leading and 
next-to-leading order, have been presented. These functions give D^ ++M (x,Q 2 ), 
but not D^ I± (x, Q 2 ). Therefore, we cannot apply them in our calculations. 

The other problem is a choice of the QCD scale parameter A and number of 
active quark flavors n/. The HT subprocesses probe the meson w.f. over a large 
range of Q 2 , Q 2 being equal to s or —u. It is easy to find that —u m in > 4.04 GeV 2 , 
while s m in > 16 GeV 2 . For momentum scales s, —t used in (42) as arguments of as 
in the LT cross section we get 

-imin > 6 GeV 2 , s min > 16 GeV 2 . 

In other kinematic domains these scales take essentially larger values. Tak- 
ing into account these facts we find it reasonable to assign A = 0.1 GeV,rif = 5 
throughout in this section. 

Results of our numerical calculations are plotted in Figs.2 — 8. First of all, 
it is interesting to compare the resummed HT cross sections with the ones ob- 
tained in the framework of the frozen coupling approximation. In Fig. 2, the ratio 
r M = (E^ T ) res /(E^ T )° for negatively charged particles (ir~,K~) is shown. In the 
computing of (E^ T )° we have neglected the meson's w.f. dependence on the scale 
Q 2 . Let us emphasize that for the kaon we have used the frozen coupling version of 
our expression (9), but not the Bagger-Gunion formula from [6], which is incorrect 
in that case. 

As is seen from Fig.2(a), r\- ^ 1 almost for all p T , whereas r K - falls from 
r K - ~ 2.75 at p T = 2 GeV/c, y = until r K - ~ 2.13 at p T = 11 GeV/c, y = and 
from r K - ~ 2.71 at p T = 2 GeV jc, y = 0.5 till r K - ~ 1.5 at p T = 11 GeV jc, y = 
0.5. For K~ this ratio demonstrates also a sharp dependence on y at fixed y/s, pr 
(Fig.2(b)). 

In all of the following figures we have used the resummed expression for the HT 
cross section. In Fig. 3 the ratio Rm =\ I I * s depicted. For all particles 
the LT cross section difference is positive Aff > 0, since £^ ~ u p (x, —t)e^, 
while Ejj^l ~ d p (x, —t)e 2 d . The smaller quark charge and the smaller distribution 
function d p both suppress [6] . The HT cross section difference may change sign 
at small pr and become negative A^ T < 0. For example, A^I < at 2 GeV jc < 
p T < 11 GeV/c for = 25 GeV, y = and at 2 GeV/c < p T < 9 GeV/c for 
x/s = 25 GeV, y = 0.5. Only at the phase-space boundary pt > 11 GeV jc in the 
first case or at p^ > 9 GeV jc in the second one E^ 7 " > E^T. Therefore, we plot the 
absolute value of R M - The similar picture has been also found for other mesons. 

As is seen from Figs. 3(a), (6) for pion and kaon the HT contribution is com- 
parable with the LT one only at pt < 3 GeV jc. We do not find a considerable 
and stable growth of HT contributions at large values of pt for all the cross section 
differences A^ T (M = ir,K), as well as, for all S^ T . Thus, E^T/E^ is small at 
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high pt for different y/s (Fig . 4(a)) . At the same time E^+/E^ is a rising function 
of p T for p T > 4 GeV/c (y/s = 14.1 GeV) and p T > 7 GeV/c (y/s = 25 GeV). 

The cross section differences Ajff and A^f = Aj^f + A^ T as functions of p T are 
shown in Fig.5. For the pion the total cross section difference A^°* in the region of 
small pt is smaller than A^ T due to A^ T < in this region (Fig. 5(a)). But for 
kaon A^* > A^ T in the same kinematic domain. For both mesons the difference 
between A M l and A^f cross sections is small. 

The rapidity dependence of Rm at y/s = 25 GeV,pT = 3 GeV/c plotted in 
Fig. 3(c) illustrates not only the tendency of the HT contributions to be enhanced 
in the region of negative rapidity, but also reveals an interesting feature of the 
HT terms; as is seen from Fig. 3(c) the ratio Rm is an oscillating function of the 
rapidity. This property of the HT terms may have important phenomenological 
consequences. In fact, in Fig.6 we have depicted A M l and Ajff versus rapidity. 
In both cases, owing to observed property of A M T (y), in certain domains of the 
rapidity interval —2 < y < 2.105 the total cross section difference is more than A^J 
and in some ones less than A^f. In the case of the kaon photoproduction 

A£* > Ajf, for -2<y< 0.3 and 1.8 < y < 2.105, 
A£* < Ajf, for 0.3<y< 1.8. 

The properties of the HT terms found in the pion and kaon photoproduction pro- 
cesses persist also in the p-meson photoproduction. But now the HT contributions 
change the whole picture of the process arising from the ordinary LT calculations. 
Thus, as in the case of the pion photoproduction, the HT terms are enhanced rela- 
tive to the leading ones and A^ T < almost for all px- But now | A^ T | takes such 
large values that it even changes the sign of the total cross section difference. That 
is, if in accordance with the LT estimations £* + > must be valid for all p T , for 
Pt < Pt we fi 11 ^ < The value of p c T depends on the process parameters, 
as well as on the p-meson w.f. used in calculations. At p t w p T we have £*+ xs E*°i. 

Our results are shown in Fig. 7. For the parameters indicated in the figure a 
critical value of px is: p T1 ~ 5.05 GeV jc for CZ w.f., and p T2 — 6.25 GeV ' jc for BB 
w.f. In all kinematic domains the HT contributions found using BB w.f. exceed 
the ones obtained by applying CZ w.f., that is, | A^ T (BB) \>\ A^ T (CZ) \. For 
example, the ratio | A^ T (B B) / A^ T (C Z) \ equals to 2.39 at = 25 GeV, p T = 
5 GeV/c, y = 0, or to 2.63 at ^ = 25 GeV, p T = 3 GeV/c, y = -1. Our 
results confirm the conclusion made by the authors in [13] concerning a possibility 
of increasing the rate of the production of transversely polarized p-meson. Similar 
pictures persist in Fig.8, where A^ T and A*°* are depicted as functions of the 
rapidity y. In Fig.8(a), for the process parameters y/s = 25 GeV,pr = 3 GeV/C 
we have: in domain V (—1.74 < y < 1.3) the total cross section difference for CZ 
w.f. is negative, in I (—1.5 < y < 1.5) - A t ° t (BB) < 0. In Fig. 8(b) the same 
is shown for y/s = 25 GeV, p T = 5 GeV/c. In two other regions lying outside of 
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/(/') the S*°+* exceeds E*°_* (for -2 < y < -1.5, £*"*(££) and for -2 < y < -1.74, 
S*°*(CZ) are negligible and are not shown). 

It is worth noticing that in [6] the authors considered the p-meson photoproduc- 
tion at the same process's parameters and predicted E*°* < at pr < 3 GeV/c, but 
could not find similar effects for S*°* in dependence on the rapidity. Our investiga- 
tions prove that < at pt < Pt and p c T well into deep perturbative domain. 
We have also demonstrated that the same phenomenon exists for y\<y< y%. 

6 CONCLUDING REMARKS 

In this work we have calculated the single meson inclusive photoproduction via 
higher twist mechanism and obtained the expressions for the subprocess jq — > Mq 
cross section for mesons with both symmetric and non-symmetric wave functions. 
For the calculation of the cross section we have applied the running coupling con- 
stant method and revealed IR renormalon poles in the cross section expression. IR 
renormalon induced divergences have been regularized by means of the principal 
value prescription and the resummed expression (the Borel sum) for the higher 
twist cross section has been found. Phenomenological effects of the obtained results 
have been discussed. 

Summing up we can state that: 

i) for mesons with non-symmetric w.f. in the framework of the frozen coupling 
approximation the higher twist subprocess cross section cannot be normalized in 
terms of a meson electromagnetic form factor; 

ii) in the context of the running coupling constant method the HT subprocess cross 
section cannot be normalized in terms of meson's elm form factor neither for mesons 
with symmetric w.f. nor for non-symmetric ones; 

iii) the resummed HT cross section differs from that found using the frozen coupling 
approximation, in some cases, considerably; 

iv) HT contributions to the single meson photoproduction cross section have im- 
portant phenomenological consequences, specially in the case of p-meson photopro- 
duction. In this process the HT contributions wash the LT results off, changing the 
LT predictions. 
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FIGURE CAPTIONS 



Fig.l Feynman diagrams contributing to the higher twist subprocess — > Mq. 
Here p and P are the hadron h and meson M four momenta, respectively. 

Fig.2 Ratio r M = (E^ T ) re 7(E^ T )°, where (Z M T ) res and (S^ T )° are HT con- 
tributions to the photoproduction cross section calculated using the running and 
frozen coupling approximations, respectively. The ratio is depicted as a function of 
Pt a), and of the rapidity b). 

Fig. 3 Ratio Rm —\ A^/A^f | for the pion a), and for the kaon b) at fixed 
rapidity y=0. In c) Rm is plotted as a function of y for the pion (dashed curve) 
and for the kaon (solid curve). 

Fig. 4 The dependence of the ratio S^ T /E^ T for K + and K~ on a) and on 

yb). 

Fig. 5 The cross section difference is shown at fixed rapidity for pions a), 
and for kaons b). For the curves 1 the process cm. energy is y^s = 14.1 GeV, for 
the curves 2 - y/s = 25 GeV. 

Fig. 6 The cross section difference A M as a function of the rapidity for pions a); 
for kaons b). 

Fig. 7 A p for p-meson. The solid curve describes A^ T , whereas the dashed curves 
correspond to A p °*. The long-dashed curve has been obtained using the CZ w.f., 
the short-dashed one- BB w.f. In the domains I(BB w.f.) and I'(CZ w.f.) the 
absolute value of | E p °* | or E*°* - E p °+* is plotted. 

Fig. 8 A p dependence on the rapidity at ^/s = 25 GeV, pt = 3 GeV/c for a); 
at y^s = 25 GeV,p T = 5 GeV jc for b). The solid curve corresponds to A^ T , the 
long-dashed and short-dashed curves describe A*°* obtained using CZ and BB w.f., 
respectively. In regions I(BB w.f.) and I'{CZ w.f.) the cross section difference 
E*°* - E to + * is shown. 
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